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Abstract. The quantised Josephson junction equation that underpins the behaviour 
of charge qubits and other tunnel devices is usually derived through cannonical 
quantisation of the classical macroscopic Josephson relations. However, this approach 
may neglect effects due to the fact that the charge qubit consists of a superconducting 
island of finite size connected to a large superconductor. We show that the well 
known quantised Josephson equation can be derived directly and simply from a 
microscopic many-body Hamiltonian. By choosing the appropriate strong coupling 
limit we produce a highly simplified Hamiltonian that nevertheless allows us to go 
beyond the mean field limit and predict further finite-size terms in addition to the 
basic equation. 
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The Josephson effect pQ is still one of the phenomena that make superconductors 
such a fascinating area of study. For instance, after some 40 years of intensive basic and 
applied research there are still new features of the coherent tunneling of Cooper Pairs 
coming to light in connection with junctions involving small superconducting grains. 
The new complications arise from the fact that when one side of a Josephson junction 
is sufficiently small for the charging energy to be relevant, quantum interplay between 
charging and tunneling begins to appear [2]. The standard approach to describing 
these junctions is to take the classical equations of motion of the superconducting phase 
difference (po across the junctions and apply cannonical quantisation rules to (pu as 
a 'position' variable. Of course, the classical Josephson equations for such a phase 
difference are first derived from microscopic theory [1], and so this standard approach 
represents a 're-quantisation' of 'classical' equations that were in turn derived from 
quantum mechanical microscopic theory using a mean field theory that does not take 
into account the charging energy of the island. Consequently, in such approaches the 
description of quantum fiuctuations is at best semi-phenomenological. In what follows 
we examine the limitations of the above procedure on the basis of a simple model 
which permits an exact treatment of a superconducting island coupled weakly to a bulk 
superconductor. 

To motivate our interest in the problem we note that much current experimental 
and theoretical attention is focused on nanoscale superconducting grains coupled to 
large superconductors as such 'Cooper Pair Boxes' are becoming realistic candidates for 
being useful qubits in Quantum Information devices [31 HI El E] . Starting with the work 
of Nakamura et al. |7j, over the last few years there have been a number of impressive 
experiments [H [HI [TOl [HI [III [13] demonstrating appropriate charge qubit behaviour and 
macroscopic tunnelling. As experiments continue to improve, it is now pertinent to re- 
examine the standard approach to describing quantum fiuctuations in superconducting 
charge qubits and related systems. Clearly, in the course of such investigation one 
would expect to reproduce the basic quantum phenomenology from a fully microscopic 
approach within a well controlled approximation, as there is already good experimental 
support for this. Nevertheless, a generalised theory will also yield new additional terms 
due to the finite size of the superconducting islands and for future experiments they could 
have significant consequences. To shed light on these, we examined a simple microscopic 
model of a Cooper Pair box, showing how the familiar phenomenology emerges, along 
with new finite size effects. 

1. Quantising the Josephson Relations 

In the interest of clarity, we start our discussion by recalling, briefly, the usual 
phenomenological approach to the problem at hand. The standard way to obtain the 
Hamiltonian describing a small superconductor connected through Josephson tunnel 
junctions to a bulk superconductor is by starting with the Josephson equations for (pjj 
and V, the difference in phases of the two superconducting regions and the voltage 
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across the junction, 

/ =Jc sin (1) 
d4>D _ 2eV_ , . 

dt ~ h ^ ^ 

where Ic is the critical current of the junction. Although derived from a quantum 
mechanical microscopic treatment, as evidenced by the appearance of h, these equations 
can be regarded as classical equations of motion. We follow the standard procedure 
for canonical quantisation and first find the Lagrangian that leads to these equations. 
Namely, we take, 

Ih'C fd<Pn\\hIc , 

L = — — ; — H COSffln (3) 

2 4e2 \ dt ) 2e ^ ^ ' 

where we have introduced the total island capacitance, C. If we choose the phase 
0D to be the canonical position variable, we can identify the canonical momentum 

TT = c}£/(90£), 

h^C ■ hCV 
^= 4^'^^ = ^ = ^(^-^^)' 
and it is seen that the phase of the condensate and the excess number of Cooper pairs 
on the island, {N — Ug), are conjugate variables OH]. The term rig represents an applied 
gate voltage (in dimensionless units) and so the canonical momentum vr = (A^ — Ug) is 
effectively the charge on the device, viewed as a capacitor, in units of 2e. Note that if the 
total charge on the system is zero, this can be rewritten in terms of the charge difference. 
To quantise the system, we introduce the commutation relation between conjugate 
variables [(po, ^r] = ih and note that this can be satisfied by writing vr = —ih d/d(f)o 
(keeping the gate voltage explicit). Then the Hamiltonian, H = vt^d — C, is given by, 

H = Ec h-^ rigj - Ej cos 0d- (5) 

where the charging energy is given by Ec = 2e^/C and Ej is defined as Ej = hlc/^e. 
The Schrodinger's equation for the amplitude i'{(pD) is then, 

Hi;i^n)=thj^iji^D). (6) 
Evidently, the probability that the phase difference takes on a certain value is given by 

This is the desired standard quantum description of the Josephson Junction. In the 
remainder of this paper, we show how the above quantised Josephson junction equation 
can be rederived directly from the microscopic theory in a way that includes finite size 
effects. 



2. Finite Superconductors As Spins 



We wish to produce a description of a finite superconductor that is simple enough to 
solve exactly but retains properties due to its finite size. Specifically, we wish to be able 
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to capture effects that go beyond the mean field approximation. Our starting point is 
the well-known BCS Hamiltonian [T^ : 

H = Y^ ekC^ ^Ck, a-^Vk, k' cl^cl^^c-k'ick'h (7) 

k, a k, k' 

where c|. ^ and Ck, a create and annihilate electrons, respectively, with spin a in the 
state k with energy and the matrix element Vk, k' describes an attractive two body 
interaction. As we are discussing a finite superconducting island, the label fc, ] does not 
refer to a free electron wavevector but to a generic single-electron eigenstate, with — /c, | 
representing the corresponding time-reversed state. 

A common approximation to this equation is made by assuming the pairing 
potential Vk^ k' is equal for all k, k' in a region around the Fermi energy determined by 
the cutoff energy hujc and zero outside this region. That is, Vk^k' = ^ for |efc — e^l < hujc 
and Vk^ fc' = otherwise. This greatly simplifies matters whilst retaining the essential 
physics. We now adopt a similar philosophy in making a further approximation, and 
take all the single electron energy levels within the cutoff region around the Fermi 
energy to be equal to the Fermi energy ep. 

The interaction term, Vk,k' acts only within the cutoff region around the Fermi 
energy. Outside this region the Hamiltonian is diagonal and trivially solved. Writing 
the single electron energy as ek = ep + {ek — ^f), we note that within the cutoff region 
\^k — ^f\ < and thus if V ^ fkJc, then |efc — e^l <C V and we can discard the 
variation of e^. Thus in this strong coupling approximation, our Hamiltonian becomes, 

H = epJ2 4, a^k, <r-Vj2 cl^clf^^c.k'iCk'^ (8) 

where the dashes on the sums indicate that they are only taken over states within the 
cutoff region. 

Although this caricature of a realistic Hamiltonian represents an uncontrolled 
approximation, we will show that it allows us to derive a Josephson junction equation 
that goes beyond mean field, and that the results it produces agree in the strong coupling 
limit with known results in two important cases. Namely, the mean-field solution of 
this Hamiltonian agrees with the BCS solution, and the exact solution agrees with the 
Richardson solution [15]. 

It should be noted that although many superconductors can be described as 
having strong coupling, the BCS Hamiltonian is not necessarily appropriate for their 
description. Equation [H] is hence not intended as a description of this particular class 
of superconductors, but rather as a generic model that, although simplified, allows an 
exact solution and a description of the physics we are trying to capture. 

As a consequence of the above simplifications equation [8] can now be written in 
terms of the three operators |16j . 

k 
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S = J2^'klCkV (9) 
k 

Note that these operators obey the commutation relations, and therefore the algebra, of 
quantum spin operators of size 1/2, where / is the number of levels in the cutoff region. 
Thus the main result of this section is the effective Hamiltonian, 

Hsp = 2{eF - fx) (^S^ + - VS+S-, (10) 

where we have introduced a chemical potential /i to describe coupling to a reservoir. 

3. Exact solution 

The eigenstates of equation [TOl are the eigenstates of the spin operator S^, ||,mAr), 
where the component of the spin along the Z axis is given bymAr = N — 1/2 and 
denotes the number of Cooper pairs on the island. The eigenenergies corresponding to 
these eigenstates are. 

En = 2{eF ~ IJ,)N -VN{1 - N + 1). (11) 

A chemical potential allows us to specify the average number of Cooper Pairs on the 
island in equilibrium. In the case of the exact solution, where is a good quantum 
number, this means choosing a state with a particular value of A^ to be the ground state. 
We choose fi so that the ground state is the state with a chosen value of A^, which we 
label N. The eigenenergies E^ therefore become, 

En= -VN{2N - N), (12) 

and the ground state is ||,m^), with energy. 

Eg, = - VN\ (13) 

We can also easily see that although the pairing parameter (5+) = in all eigenstates, 
i.e. there is no symmetry breaking, we still have fluctuations as expected for a finite 
superconductor which are given by {S^S^) = N{1 — A^ + 1) for a general eigenstate N 
and 

(|,mjv|^+^-||,m^) =iV(/-iV + l) (14) 

for the ground state. We also see that the operator S'^ which couples eigenstates, 
corresponds (when appropriately normalised) to the quasiparticle creation operator for 
the system. 

4. Comparison to Standard Results 



To generate confidence in this simple model, we compare its solutions to the solutions of 
the full Hamiltonian (equation [7j) in two ways. First, we find the mean field solution and 
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compare it to the BCS results. Second, we compare the exact solution to the Richardson 
solution in the appropriate limit. 

4.I. The Mean Field Solution 

The mean field approximation arises from the assumption that the operators remain 
close to their expectation values, {S"^). We write S"^ = (S"^) + (S*^ — {S"^)), and discard 
terms to second order or higher in {S^ — (S^)). Writing V{S^) = A, we find the mean 
field Hamiltonian for our model, 

Hmf = 2(e^ - /.) (^S^ + - AS+ - A*S- + ^ (15) 

Apart from the constant term, this is a linear combination of the spin operators , 
and and is therefore proportional to the projection of a spin operator on an unknown 
direction specified by the unit vector n. Thus denoting S.n by S?, we may write, 

Hmf = 7^1 + ^ + 2(eF - /i)^, (16) 

and therefore the problem of diagonalising equation [15] is equivalent to the problem of 
rotating the axis of quantisation for our effective spin operators. Requiring that the 
commutation relations [S?, S^] = hold for spin operators in the frame of reference 
where the axis of quantisation is along n determines both an expression for and the 
value of 7, 

2A A A* - 

St = — [s^ - — s+ - — s- 

" 7 V 2eF-7 2eF + 7 

7 =2v/(^^)2 + |Ap (17) 

where ^f = ^f ~ and we note that 7 = 2Ef, the energy of a Cooper Pair evaluated at 
the Fermi energy. The ground state of our Hamiltonian can now be trivially found, as it 
corresponds to the m = —1/2 eigenstate of the operator S?. Recalling that a maximal 
m state of a spin operator pointing in one direction is a spin coherent state [T7j in any 
other we find, 

1 (a*g+)^||,mo) 

= /nil 11(1 + "*4Ai)li"^°) (18) 

As one might expect, the second line is a way of writing the BCS ground state 
wavefunction in the limit where all the levels have equal probability of occupation, 
i.e. Uk/vk = a* for all k. Making use of equation [16] and the fact that 5'^|a) = for 
the ground state, we find that a = {^p — Ef)/ A. To complete the calculation, we need 
to self-consistently determine the values A and /x, which is relatively simple in the spin 
model and gives, 

|Ap =V'^N{l-N) (19) 
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ep-^ = V{l/2 - A^) (20) 

HsMF\a)= -riV» (21) 

where N is the average occupation of the island. We find equations [H]- EH are exactly 
equal to the expressions found if we were to solve the full BCS equation and then take 
the weak coupling limit (this is shown in the appendix of |1H])- Comparing the mean 
field solution to the exact solution, we see that, surprisingly, the exact (equation [T3l) 
and mean field (equation [2T]) ground state energies are identical. However, we have a 
non-zero pairing parameter (5*") = A/^ and the expectation value of the mean field 
coupling term, 

lA 



a 



AS+ - A*S^ + 



V 

neglects the quantum fiuctuations present in the exact solution 



2 

a ) = -VN{l-N), (22) 



- vs+s- 



^-,mj^) = -VN{l-N+l) (23) 



2' 

in much the same way a classical spin neglects the fiuctuations present in a quantum 
spin, i.e. the eigenvalues of S"^ are S{S + 1) and not the classical values S"^. 

4-2. The Richardson Solution 

Unbeknownst to the condensed matter community for many years, there exists an exact 
solution to the BCS Hamiltonian (equation [7]) for finite superconductors, first discovered 
in 1963 by Richardson |15 1 [T9l [20] in the context of nuclear physics. It has been shown 
that this solution reproduces the BCS result in the bulk limit, but it is difficult to work 
with for any island occupied by more than a few Cooper pairs. The Richardson solution 
requires the introduction of operators that diagonalise the full (i.e. not mean field) BCS 
Hamiltonian, 

N 

H = J^Ej.BjuBj. (24) 
i/=i 

fe 

where the sum in equation [2l] runs over u up to the total number of Cooper pairs on 
the island. The parameters Ej^ are found by solving the equations, 

1+ =Vy , (26) 

Ejr, — Ejy ^ 2efc — Ejy 

for all V. 

Whilst the usual BCS theory has an essential singularity at ^ = 0, the theory is 
well behaved near XjV ~ 0. Thus, following Altshuler et a/.,j21j we expand equation [261 
in powers of hjJcjV . Using ~ ep leads to, 

1 \ ^ 2 \ — \ 2(e;j — ep) /otn 

+ ^ Ej, - Ej, ^ {Ejr, - lep) + ^ (Ej, - lepY ' ^ 
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Now discard the second term on the right as neghgible, multiply by Ej^j — 2eF, and sum 
over the parameters Ej^^ to obtain, 



where the double sums over rj, v have either vanished, or gone to N{N — 1) due to 
symmetry. Finally, we recall that the energy of the Richardson solution is given by 
a sum over and rewrite equation [28] to get the energy of an island containing N 
Cooper pairs: 



As heralded in the introduction this result matches the exact energy of the spin 
Hamiltonian as given in equation [H] for /i = 0. 

Thus we have shown that although we have made a significant approximation to the 
Hamiltonian, the results thereby derived are consistent with results obtained by solving 
the full system in either the mean field approximation or exactly, and then taking the 
appropriate limit. 

5. Phase Representation of the Spin Operators and S* 

The preceding sections have established our model of a finite superconducting system 
as a large spin, as given in equation [TOl We shall now go on to show how this model can 
be used to derive a phase-representation description of the Josephson effect in a system 
comprising a small island coupled to a larger piece of bulk superconductor. 

We wish to convert to a representation in terms of the continuous phase variable </>, 
i.e. convert from ket notation to wavefunction and differential operator (such as 
^) notation. Thus, a ket will become a wavefunction and the differential 

operator must be consistent with this. Defining the state |0) = (27r)~^/^ S rriAr), 
we find that the wavefunction corresponding to ||,m7v) is (27r)^-'^/^e~*'^^. We can then 
examine how the operators act on this wavefunction. 




(28) 



E^ = 2epN - VN{1 -N+l). 



(29) 



S^(0||,mjv) 



(0|(Ar-l)|i,mAr) 
(Ar-l)(0|l,m^) 




p-i' 




(30) 



Similarly, we find for the raising operator. 



5+(0||,m^) = {ct)W[N + l){l-N)\lmM) 





(31) 
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S-{(P\^,mr,) = (0|v/iV)(/-iV + l)|i,m^) 

S-^{(j)) = e'\ i-^ 6 - ^T?r + l] ^(</>) (32) 



Collecting the differential forms for the operators and rewriting and S into a more 
convenient form leaves us with, 

S^- N— - - 
^ ~ ^'J6 2 




(33) 

The form of the raising and lowering operators can also be derived by requiring that the 
commutation relations for quantum spins are enforced. We see that it is the — |^ 
terms in the S~^, S~ operators that take into account the finite size effects and ensure 
that [S+,S-] ^ 0. 

Writing the operators in this form allows us to take the large size (/ —>■ oo) limit. 
In taking this limit we assume that 5^ ^ |. In the superconducting language, this 
corresponds to only states close to half filling being occupied. Specifically, we assume, 

L,m^\^lj)\r^Ofor\N-i\>(^-y, (34) 



a condition which is fulfilled for coherent states with set close to |. When this is 
true, we can expand the square root in (?^-|)//. We see that the leading order terms 
give 5*^ = le^**^, and we regain the semiclassical large-size limit for which [S^, S^] = 
as discussed by Lee and Scully [16]. 

6. Quantised Josephson Junction Equation 

We can now use the forms of the operators given in equations [33] to write down the 
quantised Josephson equation. We begin with a Hamiltonian that describes a finite 
superconducting island coupled to a superconducting reservoir, 

H = Hi + Hr + Hc + Ht (35) 

where Hj and Hr are the BCS Hamiltonians on the island and reservoir respectively, 
and we introduce Hamiltonians representing the charging energy of the island, 

Hc = ^{Nj-n,Y (36) 
and the tunnelling between island and reservoir, 

Ht= -TY^ c\c\c-qCq + c\c^_qC-kCk- (37) 

k,q 
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where T is the standard tunnelhng matrix element for Cooper Pairs [22], which we 
assume for simphcity to be real and independent of k, q. If we make our strong-coupling 
approximation and assume that all the electronic energy levels can be considered equal, 
we can write these Hamiltonians in terms of spin operators, as follows: 

Hr = 2{eFR — I^r)S^ — VrS^Sj^ 



2 



Hc = ^{Sf + lj/2~n,r 

Ht= - T (S+S^ + . (38) 

Inserting these expressions into equation [35|, we obtain 

H = E'aiSf - - T{SiS^ + SJS^), (39) 

where we have incorporated the terms from Hj linear and quadratic in Sf into the 
renormalised charging energy and gate voltage represented by E'^j and n'g respectively. 
In the limit that both the reservoir and the island can be considered infinite, we regain 
the standard form for the quantised Josephson junction Hamiltonian, 

d ,V IrIi 



(40) 

suggesting that equation [6] can be considered as a large-size limit where the finite size 
of the island can be neglected. However, we are now able to obtain, using equation 
the next terms in the series expansion, 

# = E^(^^A_^;^ _t||(/, + 1)cos(0,-0^) 

2 f. d li^ ^ 



We find that the new terms involve products of both the phase and the charge operators. 
Thus the Josephson tunnelling term effectively depends on the island charge. Making 
an analogy with a particle in a potential, with a position corresponding to the phase, 
we see that the extra terms in equation Hi] can be thought of as a velocity-dependent 
potential. This effect also breaks the periodicity of the island energy with n'^, i.e. the 
energy now depends on the absolute value of n^, rather than merely its value modulo 
1. In this derivation we have assumed that the single electron energies are all equal, 
and thus the occupations u^/vk are all equal. However, these occupations maintain a 
similar order of magnitude when different, and so we would expect the extra terms in 
eq. HI] to be of a similar size when the strong-coupling approximation is relaxed. 

We can make an estimate of the size of these effects for an island of a given size by 
calculating Ij, the number of electrons within the cutoff region, by comparing the level 
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spacing to the Debye energy. For an island with a volume equal to that described in [7j , 
we find that // ^ 6 x 10^, and thus the additional terms in equation UT] (proportional 
to l/lj are unlikely to be significant. If we consider instead a nanograin of the type 
described in [191 ISOl [23], we find that // ~ 400, and thus the extra terms may be 
relevant. 

7. Conclusions 

We have shown how the quantum Josephson Junction equation, usually derived by re- 
quantising the mean field equations of motion, can be directly derived from a microscopic 
description of a superconducting island. We used a simplified Hamiltonian in which the 
energy of the individual microscopic electron levels is considered equal that allowed 
an exact solution to be found. We have shown how a mean field approximation leads 
to a solution that corresponds to a spin coherent state, which is the BCS state in 
the appropriate limit. As well as illustrating how the familiar phenomenology emerges 
through the mean field approximation, we showed we can describe effects beyond the 
mean field, such as quantum fiuctuations. We went on to rederive the Josephson 
Junction equation and describe size dependent corrections to the familiar terms. 
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